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Most quantum computation schemes propose encoding 
qubits using natural two level systems such as spin 1/2 parti- 
cles. Others exploit only two states of a larger discrete Hilbert 
space such as the energy excitation levels of an ion. In the 
paper "Encoding A Qubit In An Oscillator" Ql, Gottesman, 
Kitaev, and Preskill (GKP) described an alternative method 
for encoding a qubit in the continuous position and momen- 
tum degrees of freedom of an oscillator. Because the qubit 
is encoded in an infinite dimensional Hilbert space, a natural 
error correcting scheme arises. It can correct errors due to 
shifts in the oscillator's position or momentum (i.e., small ap- 
plication of the displacement operator). When the oscillator is 
chosen to be a single mode of the electromagnetic field, fault 
tolerant computation can be performed by use of only phase 
shifters, beam splitters, squeezing, photon counting, and ho- 
modyne measurements. 010]. However, state preparation 
requires nonlinear interactions. 

The GKP scheme constitutes a type of linear optical quan- 
tum computer. Other schemes for such computers are based 
on the proposals of Knill, Laflamme, and Milburn (KLM) 
H |1 13 and of Ralph et al. @]. Of these the KLM pro- 
posal has received most of the analysis H.I7II8II9I [Toll . How- 
ever, which (if any) of these three schemes is most suitable 
for large scale quantum computation is still unknown. Thus, 
further analysis of all three schemes is required. This paper 
is devoted to an analysis of the GKP quantum computer. We 
consider errors only in the input qubit and ancilla states, and 
we obtain a threshold for the maximum size of position and 
momentum shifts that the error correcting scheme can always 
correct. We then use this threshold to constrain the quality of 
the input states and estimate the number of photons the input 
states must contain. 

We first give a brief review of qubits in the GKP scheme. 
One may use any type of oscillator to represent a qubit, but 
for the purposes of this paper we choose the oscillators to be 
single modes of the electromagnetic field. Let di be the pho- 
ton annihilation operator of mode i, d\ the creation operator, 



Xi — 1/ ' ^/2(d\+di) the x-quadrature, and pi — i/\/2{d\—di) 
the p-quadrature. Let \x) denote eigenstates of x, and \p) 
eigenstates of p. We represent the logical qubit as 
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and the logical 1 qubit is 
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In the x-quadrature these states are infinitely long combs of 
delta functions, with the | 1l) displaced a distance y/ir from 
the |0l) state. Such states are clearly unphysical, having an 
infinite energy expectation value, and they are not normaliz- 
able. 

This encoding can protect the qubits from errors due to dis- 
placements in the x- and p-quadratures, described by the op- 
erators e _m P (shifting the x-quadrature a distance of u) and 
e -zvx (shifting me p-quadrature a distance of v). These shift 
operators form a complete basis, so any error superoperator £ 
acting on the density operator p of a single oscillator may be 
expanded as 1 1 ] 



£{p)= / du&vrlu'&v' F(u,v,u' ,v')e- luf, e 
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If the distribution F(u,v,u' ,v') is sufficiently concentrated 
near the origin, then this encoding allows us to correct £ and 
recover p. 

Measurement of the qubits can be accomplished by use of 
homodyne detection of the x-quadrature. Results in which the 
x-quadrature is measured to be closer to an even multiple of 
y/ir are registered as the detection of |0l), and measurements 
closer to an odd multiple are registered as Therefore, a 
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shift error that is larger than \pH j2 may cause an error in the 
qubit measurement. 

Because the qubit states are unphysical, we must find some 
states that can approximate the ideal qubit states and are 
physically realizable. GKP propose using states whose x- 
quadrature wave function is composed of a series of Gaussian 
peaks of width A contained in a larger Gaussian envelope of 
width 1/fc. The approximations of the and 1 logical states 
are 



(z|6> - N 



oo 

E 

s— — oo 



and 



(6) 



(a) 



l(x|0)| 2 



1 

0.8 
0.6 
0.4 
0.2 



-4^ -2^ 



11 



xl) 



N1 £ 



e -!((2s+l)fcv^) 2 e 



*0 



(7) 



where No and Ni are normalization factors. In the limit that k 
and A are both small, No ~ iVj = \/2fc. We show examples 
of these approximate qubit states in Fig. First, note that 
these states are not orthogonal, and there will be some proba- 
bility of mistaking a state for a 1 state (and a 1 for a 0). This 
probability is equal to the probability that the ^-quadrature of 
the state is measured to be closer to an odd multiple of y/n. 



Po^i= E / _ dx\(x\0)\ 
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In the limit of small A and k, this expression can be approxi- 
mated as 
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For A = fc = 0.5, P ^i ~ 0.01, and for A = k = 0.25, 
P ^i ~ 10~ 6 . For small A and k, P x ^o ~ Po->i- How- 
ever, we desire not only that measurements can distinguish 
between the and 1 states, but also that an error correcting 
procedure can reliably correct errors. This will place stricter 
requirements on the approximate qubit states. 

We now examine the error correcting circuits. We imagine 
that the qubit \Ql) initially exists in a superposition of the 
|0l) and states in mode 1. The qubit then receives errors 
resulting in shifts of u\ in x\ and v\ in p\. We begin the cor- 
rection procedure by repairing the shifts in x\. This requires 
an ancilla qubit prepared in the state |+ L ) = |0l) + in 
mode 2. Of course, the ancilla may also be subject to errors, 
resulting in shifts u 2 and v 2 - After the errors, the state of the 
system is 
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To correct the qubit's errors, modes 1 and 2 are sent through 
the network pictured in Fig. [2] The two damaged modes meet 



FIG. 1: Example approximate qubit states, (a) shows the x- 
quadrature wave function of 1 0) , and (b) shows j 1) . For each of these 
we have chosen A = k = 1/4 



in a beam splitter, which performs the transformation 
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in the x-quadrature wave function of the two modes. After 
mode 1 exits the beam splitter, we apply the squeezing oper- 
ator S(y/2), where S is defined by S(q)\x) = y/q\x/q). We 
now measure the x-quadrature of mode 2 and obtain the re- 
sult X2 — 1/V2 (ny/n — u\ — U2) , where n may be any inte- 
ger. This measurement provides some information about the 
errors shifting the x-quadrature of the qubit. The errors can 
be (partially) corrected by applying the displacement operator 
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Here the Mod function has the range [— s/n, s/n). 

This entire sequence of operations - errors and correction 
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FIG. 2: Procedure to correct x shifts acting on the qubit. The circuit 
shows errors acting on the qubit and ancilla. Then the qubit and 
ancilla modes meet in a beam splitter, mode 1 is squeezed, the x- 
quadrature of mode 2 is measured, and a correcting shift of s(x2) is 
applied to the qubit. 



qubit = QJ 
ancilla = , ) 
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FIG. 3: Procedure to correct p shifts acting on the qubit. The circuit 
shows errors acting on the qubit and ancilla. Then the qubit and 
ancilla modes meet in a beam splitter, mode 1 is squeezed, the p- 
quadrature of mode 2 is measured, and a correcting shift of s(p2) is 
applied to the qubit. 



procedure - transforms the qubit to the state 



e itp(u 1 ,v 1 ,u 2 ,v 2 ,n) e -i(v 1 -v 2 )xi 



xe -i('«i-5 Mod 2v^( 2 «i+ 2 «2))pi |Q L ^ (15) 

where e'^Oi^i^,^,") j s a ph ase factor independent of the 
input qubit. If |ui + w 2 | < \p^l% then the qubit is "corrected" 
to the state 



|Q L ^ > gi0(t»i,Di,«2,f2,n) e -i(fl-t>2)£l 



xe 



IQl 



(16) 



Notice that both the ancilla and the qubit's p shift errors v± and 
v 2 both appear now acting on the qubit. Also, the qubit has 
been affected by the ancilla's x shift error u 2 . In the case when 
\u\ + u 2 \ > the error correcting procedure applies the 

Pauli X operator to the qubit state, producing 



|Q L ^ j. e i4>(ui,v 1 ,u 2 ,v 2 ,n) e -i(v 1 -v 2 )x 1 
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Correcting the X error would require the use of a standard 
quantum error correcting code (concatenated with this contin- 
uous variable code). See 1 1 1] for an introduction to quantum 
error-correcting codes. 

The procedure for correcting shifts to the qubit's p- 
quadrature wave function is similar to that used for x shifts. It 
is pictured in Fig. [5] In this case we prepare the ancilla qubit 
in the state |0l), we measure the p-quadrature basis of mode 
2, we apply the squeezing operator 5(1/ V2), and the correct- 
ing shift is e -ls (P2)zi tj^ ell - ors anc [ correcting procedure 
produces the new state 
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where the phase factor e l6 ( u i, v i> u 2> v 2> n ) i s independent of the 
initial qubit state. When \v\ + v 2 \ < \f^j2 the correcting 
procedure is successful, resulting in the state 



\Q L j j, e i8(ui,vi,u 2 ,v 2 ,n) e ~ i(ui-u 2 )pi 

xe- W2£l |Q L >. 
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When vi and v 2 are too large, the correcting circuit creates a 
Pauli Z error in the logical qubit basis. 

Armed with a clear understanding of the error correction 
procedure, we can formulate a bound on the maximum error 
that this scheme can always correct. In the following we imag- 
ine that the only error source is in the preparation of the qubit 
and ancilla states. We assume that all of the operations of the 
error correcting procedure are error free. In the above analy- 
sis we have shown how errors in ancillas are transferred to the 
qubit. We would like to find a bound on the size of shift er- 
rors affecting the qubit and ancillas that ensures that the qubit 
never receives an X or Z error and that the sizes of x and 
p shifts do not grow after multiple steps of error correction. 
To accomplish this we follow the qubit's x and p shifts as it 
passes through multiple steps of correction. 

Before the error correction procedure, the qubit has errors 
iii and v\, and the first ancilla has errors u 2 and v 2 . If we first 
correct errors of the x-quadrature, the qubit will have errors 
u 2 (in x) and v\ — v 2 (in p), provided that \u\ + U2I < \f^l^- 
This result comes directly from Eq. d!6i . Now, we should 
correct errors of the p-quadrature using a new ancilla, which 
may have errors 113 and V3. The p error correction produces a 
new qubit with errors u 2 — U3 (in x) and v% (in p), provided 
that — v 2 + W3I < At this point, the first stage of 

error correction is complete. The qubit's original errors have 
been entirely eliminated and replaced by errors introduced by 
the ancillas. A second correction of the x-quadrature using 
a third ancilla with errors 114 and V4 is successful if \u 2 — 
u 3 + u 4 \ < \fn/2. The qubit now has errors U4 and v$ — 
v 4 . In this and each subsequent stage, we find that at any 
time, the qubit has errors inherited from one ancilla in the 
quadrature that was just corrected and errors from two ancillas 
in the quadrature, which should be corrected next. During the 
next correction errors from a new ancilla are introduced. The 
correction succeeds if the magnitudes of the now three errors 
is less than Therefore, we can be certain that repeated 

error correction steps will be successful if the magnitude of 
all errors shifts are smaller than t/k/Q. 

We would now like to calculate the probability that an ap- 
proximate qubit state such as |0) or |1) has shifts smaller than 
y/Tr/6 in both x- and p-quadratures. To find this probability 
we decompose the approximate qubit states in a basis defined 
by the states that are x- and p-shifts of |0l), which we may 
express as 

|tt,») = 7r- 1 - /4 e - i ^e- i ^|Q L ) (20) 
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= 7r_V4 £ e ' m2sV ^\ x = 2sy/^+u). (21) 



xe -i(v 1 -±Mod 2 ^(2v 1 +2v 2 ))x 1 \Q L ^ (18) We can express any state |t/>) of a single oscillator in this basis 
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FIG. 4: P(u,v) the probability density that the |0) approximate 
qubit state contains errors u (in its z-quadrature) and v (in its p- 
quadrature). Here we have used A = k = 1/4. 




FIG. 5: P(u,v) the probability density that the |1) approximate 
qubit state contains errors u (in its x-quadrature) and v (in its p- 
quadrature). Here we have used A = k = 1/4. 



using a wave function f(u, v) = (u, v\ip): 



du / dvf(u, v)\u, v). (22) 



The limits of the above integral are chosen to match the pe- 
riodicity in |0l). The probability density for a state having 
shifts u and v from |0l) is simply 



P{u,v) = \f{u,vf 



(23) 



Examples of P(u, v) are shown in figures0]and[5] 

The probability that the |0) state has errors less than 
(which guarantees that its errors may always be corrected) is 

6 du I dvP(u,v). (24) 

6 6 

We plot P no en-or in Fig. |6] Achieving P no enor — 0.9 requires 
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FIG. 6: Here we plot P no error, the probability that an approximate 
qubit state has shifts in both x- and p-quadratures less than \/tv/6 as 
a function of A = k. This shows the case of |0), but the |1) plot is 
indistinguishable. 
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FIG. 7: Here we plot the mean number of photons contained in an 
approximate qubit state as a function of the probability that such a 
state has a shift error larger than \fn/ 6. 



A = k = 0.214, and P no enor = 0.99 requires A = k = 
0.149. 

We generally believe that protecting an oscillator from de- 
coherence usually becomes more difficult when the oscillator 
contains large numbers of photons, so we are motivated to cal- 
culate the mean number of photons in the approximate qubit 
states. A crude estimate of this quantity for small A and k is 
given by 
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In Fig. Qwe plot the mean number of photons contained in an 
approximate qubit state as a function of P e rror = 1 — Pno error 
using an exact expression. We find that an approximate qubit 
state with P no en-or = 0.9 must contain (n) = 10.4 photons, 
and to achieve P no error = 0.99 requires (n) = 22.1 photons. 

Production of these states is likely to be very challenging. 
We are aware of four proposals for production of qubit states 
for the GKP quantum computer. The first is from the original 
work of GKP 1 1], in which they propose the use of a two- 
mode Hamiltonian of the form xifi2- This might be achieved 
by coupling an optical mode with a mirror that may exist in a 



5 



coherent superposition of position states. For a discussion of 
the possibility of such experiments see fl2lfT3ll and their refer- 
ences. The second proposal is from Travaglione and Milburn 
0. They describe a method that prepares the qubit states in 
the oscillatory motion of a trapped ion rather than the photons 
in an optical mode. The third proposal is by Pirandola et al. 
llLHl and discusses the preparation of optical GKP states using 
a two mode Kerr interaction (described by a Hamiltonian of 
the form h\fi2) followed by a homodyne measurement of one 
of the modes. The same authors also describe a fourth method 
for GKP state production in lU6lfT7ll . This method would pre- 
pare the GKP state in the motion of a neutral atom that in- 
teracts with a photon field confined to a high finesse optical 
cavity. They show that this may be done so that the atom is 
trapped in the cavity, or the atom may be free to pass through 
the cavity. State preparation is currently the most problematic 
aspect of the GKP quantum computation scheme, and it is an 
area in need of more experimental effort. 

In this paper we have investigated error correction in the 
GKP quantum computer, considering errors in both qubit and 
ancilla states. Provided that all input states have displacement 
errors in their x- and p-quadratures less than t/tt/6, the er- 
ror correction schemes will operate successfully. Because it 
is impossible to prepare perfectly error free GKP states, we 
have calculated the probability that approximate GKP states 



contain an error larger than y/n/ 6. An approximate GKP state 
with error probability less than 0.1 must contain a mean num- 
ber of photons greater than 10.4. 

There are some cases in which states with errors larger than 
will not cause logical qubit errors. For example, the 
qubit's and the ancilla's errors may actually cancel one an- 
other. Because of effects like this, our y/n/ 6 bound may be 
refined with more detailed analysis. We have not considered 
full fault tolerant computation in the presence of noisy logic 
operations, nor have we considered the effects of phase er- 
rors or photon absorption. Phase errors and photon absorption 
are likely to be primary error sources for this (or any) optical 
quantum computer. Because linear superpositions of displace- 
ments span the space of possible single qubit errors, one may 
correct phase and absorption errors using the GKP error cor- 
rection scheme. However, phase errors affecting states with 
large numbers of photons correspond to large displacements. 
This fact highlights two competing influences: (1) To lower 
the number of intrinsic errors in each approximate qubit state, 
we must prepare states with large numbers of photons. (2) To 
reduce the displacements caused by phase errors, we should 
prepare states with fewer photons. A detailed analysis of these 
constraints requires further study. 
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